Tickling the CMB damping tail: scrutinizing the tension between the ACT and SPT 

experiments. 
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The Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT) have recently 
provided new, very precise measurements of the cosmic microwave background (CMB) anisotropy 
damping tail. The values of the cosmological parameters inferred from these measurements, while 
broadly consistent with the expectations of the standard cosmological model, are providing interest- 
ing possible indications for new physics that are definitely worth of investigation. The ACT results, 
while compatible with the standard expectation of three neutrino families, indicate a level of CMB 
lensing, parametrized by the lensing amplitude parameter Al, that is about 70% higher than ex- 
pected. If not a systematic, an anomalous lensing amplitude could be produced by modifications 
of general relativity or coupled dark energy. Vice-versa, the SPT experiment, while compatible 
with a standard level of CMB lensing, prefers an excess of dark radiation, parametrized by the 
effective number of relativistic degrees of freedom N e g. Here we perform a new analysis of these 
experiments allowing simultaneous variations in both these non-standard parameters. We also com- 
bine these experiments, for the first time in the literature, with the recent WMAP9 data, one at 
a time. Including the Hubble Space Telescope (HST) prior on the Hubble constant and informa- 
tion from baryon acoustic oscillations (BAO) surveys provides the following constraints from ACT: 
N cB = 3.23 ± 0.47, A L = 1.65 ± 0.33 at 68% c.L, while for SPT we have N eS = 3.76 ± 0.34, 
Al = 0.81 ±0.12 at 68% c.L. In particular, the Al estimates from the two experiments, even when 
a variation in 7V c ff is allowed, are in tension at more than 95% c.L. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 
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I. INTRODUCTION 

The new measurements of the Cosmic Microwave 
Background (CMB) anisotropics provided by the Ata- 
cama Cosmology Telescope (ACT) [T| and the South 
Pole Telescope (SPT) @ have both provided new and 
exquisitely precise observations of the CMB damping tail. 

This angular region of the CMB angular spectra, cor- 
responding to the multipole range going from £ ~ 700 up 
to I ~ 3000, plays a key role in the determination of cru- 
cial parameters like the relativistic number of degrees of 
freedom N e ff, the primordial Helium abundance Y p and 
the running dn/dh\k of the scalar spectral index. 

Among those parameters, Y p can be determined un- 
ambigously assuming standard Big Bang Nucleosynthesis 
(and thus does not represent a free parameter of the the- 
ory), while dn/dhik is expected to be negligible in most 
inflationary models. On the other hand, the effective 
number of relativistic degrees of freedom N e g practically 
parametrizes the energy density of relativistic particles 
in the early Universe. In the standard scenario, with 
the three active relativistic neutrino species, a value of 
N e g = 3.046 is expected Q. Deviations from this value 
due to a non-vanishing neutrino chemical potential are 
possible but bound to be small, especially in light of the 
recent evidences for a large value of the neutrino mixing 
angle #13, see e.g. Thus a detection of N e s ^ 3.046 



would point to the presence of physics beyond the stan- 
dard model of particle physics, like the existence of a yet 
unknown particle, e.g., a sterile neutrino. 

The damping tail is not only affected by those param- 
eters but also by other physical effects generally tak- 
ing place at a much later epoch, well after recombi- 
nation. These includes, for example, the extragalactic 
foreground emission of point sources, radio galaxies, the 
Sunyaev-Zel'dovich effect and similar unresolved back- 
grounds. These foregrounds can however be well identi- 
fied by their spectral and angular dependence and have 
in general a minimal correlation with the cosmological 
parameters. 

More importantly, the CMB damping tail is affected 
by the lensing of CMB photons by dark matter clumps 
along the line of sight. This effect is linear, can be com- 
puted precisely and depends on the same cosmological 
parameters that affect the primary CMB spectrum. How- 
ever, the lensing amplitude is strictly dependent from the 
growth of perturbations. This quantity can be significant 
different if, for example, general relativity is not the cor- 
rect theory to describe gravity at the very large scales. If 
the accelerated expansion of our universe is indeed pro- 
vided not by a dark energy component but by modified 
gravity, the perturbation growth could be dramatically 
different and change the expectations of lensing (see, for 
example @ and references therein). In order to test the 
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correct amplitude of the lensing signal, one can introduce 
a calibration parameter Al , as in [7] , that scales the lens- 
ing potential in such a way that Al = corresponds to 
the complete absence of lensing, while Al — 1 is the ex- 
pected lensed result assuming general relativity. A robust 
detection of Al being different from unity would hint to 
the fact that general relativity is not the correct theory 
to describe gravity at the cosmological scales. 

The new ACT and SPT data, while broadly consistent 
with the expectations of the standard ACDM scenario, 
are indeed providing interesting hints for deviations from 
the simplest ACDM model when combined with the re- 
sults from 7 years of observations from the WMAP satel- 
lite (WMAP7, [13]). 

The SPT experiment, for example, is confirming an 
indication for a value for N e g > 3.046. This indica- 
tion, already present in the previous data release (see 
e.g. [|, @ and 0), is marginal when considering only 
the WMAP7+SPT data with N eB = 3.62 ± 0.48 at 68% 
c.L. However, it is more significant when the SPT data 
is combined with the measurement of the Hubble con- 
stant H = 73.8 ± 2.4km s _1 Mpc" 1 from the Hubble 
Space Telescope (HST) [ll[ and with information from 
Baryonic Acoustic Oscillation (BAO) data (see Table 4 
in Q), yielding a final value of N e ft = 3.71 ± 0.35. 

At the same time, the ACT collaboration presented 
a similar analysis obtaining different results. In particu- 
lar, the WMAP7+ACT data alone constrain the neutrino 
number to be N e g = 2.78 ± 0.55, i.e. perfectly consis- 
tent with the standard three-neutrino framework. When 
the ACT data is combined with HST and BAO data the 
value is higher, N c s — 3.52 ± 0.39, but still consistent 
with three neutrinos families (see Table III in [1]). 

Interestingly, this is not the only tension between the 
two datasets. If we now consider the results on the 
lensing amplitude parameter, the SPT dataset is fully 
compatible with the standard expectation, with Al = 
0.86±°;^ at 68% c.l. (see [lj]), while the ACT data sug- 
gest a 2cr deviation from the standard expectation, with 
A L = 1.70 ±0.38 at 68% c.l.. 

In this brief paper we further investigate these discrep- 
ancies by improving these analyses in two ways. First of 
all, we perform our analyses allowing both N e g and Al 
parameters to vary at the same time. As we will see, this 
let to better identify the tension between the two exper- 
iments. Secondly, we add the recent dataset from nine 
years of observations coming from the WMAP satellite 
as in [l4|. Both ACT and SPT teams used the previous 
7-year WMAP dataset in their papers and some, albeit 
small, differences are present when the updated dataset 
is considered. 

Our paper is simply organized as follows: in the next 
section we describe the analysis method, in Section III 
we present our results and in Section IV we derive our 
conclusions. 



II. DATA ANALYSIS METHOD 

Our analysis is based on a modified version of the pub- 
lic CosmoMC [H| Monte Carlo Markov Chain code. We 
consider the following CMB data: WMAP9 Q, SPT 0, 
ACT Q including measurements up to a maximum mul- 
tipole number of Z max = 3750. For all these experiments 
we make use of the publicly available codes and data. For 
the ACT experiment we use the " lite" version of the like- 
lihood [22]. Since ACT and SPT dataset are providing 
different results on the parameters, we will consider them 
separately. Thus our basic CMB-only datasets consist of 
the WMAP9+ACT and WMAP9+SPT data. 

We also consider the effect of including additional 
dataset to the basic datasets just described. Consistently 
with the measurements of the HST [ll|, we consider 
a gaussian prior on the Hubble constant Ho = 73.8 ± 
2.4kms~ 1 Mpc -1 . We also include information from 
measurements of baryonic acoustic oscillations (BAO) 
from galaxy surveys. Here we follow the approach pre- 
sented in [l4j combi ning four datasets: 6dFGRS from 
0, SDSS-DR7 from [lj, SDSS-DR9 from QJ] and Wig- 
gleZ from Ql]. 

We sample the standard six-dimensional set of cos- 
mological parameters, adopting flat priors on them: the 
baryon and cold dark matter densities Q^h 2 and fl c h 2 , 
the ratio of the sound horizon to the angular diameter 
distance at decoupling 9, the optical depth to rcionization 
r, the scalar spectral index n s , the overall normalization 
of the spectrum A s at k = 0.002 Mpc -1 . 

Since the ACT and SPT data are showing indications 
for deviations from their standard values, we also con- 
sider variations in the effective number of relativistic de- 
grees of freedom N e g and in the lensing amplitude pa- 
rameter Al as defined in Q, that simply rescales the 
lensing potential: 

C+* -> A L Cf* (1) 

where Cf is the power spectrum of the lensing field. We 
take flat priors on all the parameters; in particular, we 
take 1 < Acs < 10 and < A L < 4. 

In our basic runs, we do not consider the effect of mas- 
sive neutrinos. We perform additional runs in which we 
allow for a non- vanishing neutrino mass, parametrized 
by means of the neutrino fraction /„ = Q, u /Q c . We al- 
ways assume standard Big Bang Nucleosynthesis, so that 
the Helium abundance Y p is uniquely determined by the 
values of il^h 2 and N e s . 

Finally, in order to assess the convergence of our 
MCMC chains, we compute the Gelman and Rubin R — 1 
parameter demanding that R — 1 < 0.03. 

III. RESULTS 

As stated in the previous section, we consider the 
ACT and SPT datasets separately. We therefore 
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FIG. 1. Constraints in the Al - N c fi plane from a CMB only analysis (left panel) and including the HST prior and BAO (right 
panel). The blue contour includes the ACT data while the red contour refers to the SPT data. The line at Al = 1 indicates 
the standard expectations based on General Relativity. The line at N e g = 3.046 indicates the prediction from the standard 
model with three neutrino flavours. 




FIG. 2. Constraints in the Q.bh 2 - n a plane (left panel) and on the Q c h 2 -Ho plane (right panel) from ACT (blue contours) and 
SPT (orange contours) including WMAP9, HST and BAO data. The ACT-SPT tension is less pronounced for these parameters. 



perform the following four analyses: WMAP9+ACT, 
WMAP9+ACT+HST+BAO, WMAP9+SPT, and 
WMAP9+SPT+HST+BAO. 

In Table I we report the constraints on the considered 
parameters from each run. As we can see, the ACT and 
SPT are providing significantly different constraints on 
the N e s and A L parameters. 

In order to further investigate this discrepancy, we plot 
in figure Q] the 2-D constraints on the N c g vs A L plane 
for the CMB only case and for the CMB+HST+BAO 
analysis. 

As we can see, a tension is clearly present since the cen- 
tral values for iV off and A L obtained from WMAP9+ACT 
analysis are outside the 95% confidence level of the 
WMAP9+SPT and vice-versa. Namely, the ACT dataset 
is pointing towards a value of N c g fully consistent with 
the standard scenario of N e g = 3.046, while (as it can 



be seen from Table 1 and Figure I), preferring at the 
same time an exotic high value for the lensing potential, 
with A L larger than unity at more than 95% c.l. when 
the BAO and HST datasets are included. Considering 
the 95% confidence levels we found A L — 1.64±°;^ for 
the WMAP+ACT analysis and A L = 1.65±g-|l for the 
WMAP+ACT+BAO+HST. 

The situation is opposite for the SPT data: while SPT 
is fully consistent with Ai, = 1, N c g is constrained to a 
larger value than the standard expectation. When also 
the HST and BAO data are included, we see that not only 
a value of N c g > 3.04 is suggested at more than 95% c.l., 
but also a value of A L smaller than one is suggested at 
about 68% c.l.. 

In particular, we found that A L < 1.07 at 95% c.l. 
from WMAP9+SPT+BAO+HST while A L > 1.13 at 
95% c.l. from WMAP9+ACT+BAO+HST, i.e. for the 
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Parameters 


SPT+WMAP9 


ACT±WMAP9 


SPT±WMAP9±HST±BAO 


ACT±WMAP9±HST±BAO 


tt c h' 2 
6 

T 

n 3 

AW 

A L 

#o[km/s/Mpc] 
log(10 10 A) 

Age/Gyr 

r\SZ 

-^3000 

n CL 

-^3000 

n PS 

^3000 

A sz 


0.02264 ± 0.00051 
0.1232 ±0.0080 
1.0415 ± 0.0012 
0.088 ±0.014 
0.982 ±0.018 
3.72 ± 0.46 
0.85 ±0.13 
74.6 ± 3.7 
3.169 ±0.048 
0.736 ± 0.023 
0.264 ± 0.023 
13.14 ± 0.43 
5.8 ±2.4 
5.2 ±2.1 
19.6 ± 2.5 


0.02283 ± 0.00052 
0.110 ±0.010 
1.0412 ± 0.0025 
0.090 ±0.014 
0.969 ±0.019 
2.85 ±0.56 
1.64 ±0.36 
69.9 ±3.7 
3.174 ±0.045 
0.728 ± 0.025 
0.272 ± 0.025 
13.90 ±0.55 

1.00 ±0.57 


0.02255 ± 0.00036 
0.1274 ± 0.0075 
1.0413 ± 0.0012 
0.085 ±0.013 
0.979 ±0.011 
3.76 ± 0.34 
0.81 ±0.12 
73.3 ± 1.8 
3.185 ±0.034 
0.721 ±0.015 
0.279 ±0.015 
13.12 ±0.28 
5.8 ±2.5 
5.3 ±2.1 
19.5 ±2.5 


0.02294 ± 0.00042 
0.1178 ±0.0094 
1.0400 ± 0.0024 
0.090 ± 0.014 
0.978 ±0.014 
3.23 ±0.47 
1.65 ±0.33 
71.1 ± 2.4 
3.174 ±0.037 
0.721 ±0.017 
0.279 ±0.017 
13.55 ±0.41 

0.91 ±0.57 


Xmin 


3806.25 


3798.83 


3808.06 


3800.59 



TABLE I. Cosmological parameter values and 68% confidence level errors. The SPT and ACT datasets produce different values 
for some of the parameters, most notably N e fi and Al- 



Parameters 


SPT 


SPT 


ACT 


ACT 




±WMAP9±HST±BAO 


±WMAP9±HST±BAO 


±WMAP9±HST±BAO 


±WMAP9±HST±BAO 




0.02278 ± 0.00036 


0.02274 ± 0.00037 


0.02289 ± 0.00042 


0.02303 ± 0.00043 


n c h 2 


0.1305 ±0.0084 


0.1310 ±0.0082 


0.1163 ±0.0087 


0.1193 ± 0.0097 


e 


1.0412 ±0.0011 


1.0412 ±0.0011 


1.0402 ± 0.0023 


1.0402 ± 0.0023 


T 


0.089 ± 0.013 


0.089 ± 0.013 


0.094 ±0.014 


0.091 ±0.014 


n a 


0.989 ±0.012 


0.987 ±0.012 


0.976 ±0.014 


0.981 ±0.014 




3.89 ±0.38 


3.88 ±0.37 


3.09 ±0.43 


3.28 ±0.48 


Em„ [eV] 


0.40 ±0.22 


< 0.77 (95% c.l.) 


< 0.40 (95% c.l.) 


< 0.53 (95% c.l.) 


A L 


1.00 


0.90 ±0.15 


1.00 


1.78 ±0.38 


if [km/s/Mpc] 


72.4 ± 2.0 


72.5 ± 1.9 


69.3 ± 2.4 


70.1 ±2.5 


log(10 10 A) 


3.156 ±0.034 


3.166 ±0.036 


3.177 ±0.036 


3.162 ±0.037 


n A 


0.707 ± 0.020 


0.707 ± 0.019 


0.710 ±0.019 


0.710 ±0.020 




0.293 ± 0.020 


0.293 ±0.019 


0.290 ±0.019 


0.290 ± 0.020 


Age/Gyr 


13.12 ± 0.29 


13.10 ± 0.28 


13.75 ±0.40 


13.59 ±0.42 


n sz 

^3000 


5.8 ±2.4 


6.1 ± 2.4 






n CL 
^3000 


5.3 ±2.2 


5.2 ± 2.1 






n PS 
^3000 


19.3 ± 2.4 


19.3 ± 2.5 






A S z 






0.97 ±0.57 


0.97 ±0.57 


Xmin/2 


3808.0 


3807.5 


3802.47 


3800.59 



TABLE II. Cosmological parameter values and 68% confidence level errors for the analysis that considers massive neutrinos. 
As we can see, varying Al strongly affects the constraints on the total neutrino mass. Vice- versa, allowing for a neutrino mass 
renders the SPT value for Al more compatible with the standard value while exacerbates the problem for the ACT dataset. 



lensing parameter the SPT and ACT datasets are pro- 
viding constraints that are in disagreement at more than 
95% c.l.. 

It is interesting to note that the tension between the 
ACT and SPT datasets is clearly not limited to Al or 
N c g: also the constraints on i?o, n s , ft^h 2 and Q c h 2 ap- 
pear as quite different. The discrepancy is however less 
significant since the central values are inside the 95% 
confidence level of each analysis (see Figure [5]). 

The results discussed so far are relative to the analysis 



in which all neutrinos are considered as relativistic and 
massless. Since the SPT dataset is claiming a detection 
at 95% c.l. for a neutrino mass with Em^ = 0.48 ± 0.21 
in a WMAP7±SPT±BAO±HST analysis (see 0), it is 
clearly interesting to consider also massive neutrinos. 

In table II we present the constraint on cosmological 
parameters from the WMAP9±SPT±HST±BAO and 
WMAP9±ACT±HST±BAO datasets respectively when 
variation in the neutrino masses are included in two cases: 
varying Al and fixing Al = 1. 
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FIG. 3. Posterior distribution function for the total neutrino mass parameter Em„ from a SPT+WMAP+BAO+HST analysis 
(left panel) and ACT+WMAP+BAO+HST (right panel) in the case of fixing lensing to Al = 1 and letting it to vary. As we 
can see, if we let the Al parameter to vary the small indication for a neutrino mass from the SPT analysis vanishes. At the 
same time, letting the Al parameter to vary weakens the constraints from ACT. 
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FIG. 4. Constraints in the Al vs Em„ plane for the 
SPT+WMAP+BAO+HST and ACT+WMAP+BAO+HST 
datasets. A degeneracy is present between the two param- 
eters: larger values for Al let larger neutrino masses to be 
more consistent with the data. The SPT indication for a neu- 
trino mass is driven by the low value of Al obtained in the 
neutrino massless case. 



As we can see, while the ACT datasct does not favour 
the presence of neutrino masses, the SPT dataset gives 
Era,, = 0.40 ± 0.22 at 68% c.l. in the case of A L = 1 and 
a lower limit limit Sm„ > O.OAeV at 95% c.l.. This is 
consistent with the results reported in Q considering the 
different WMAP and BAO datasets. However, when Al 
is let to vary, the evidence for a neutrino mass vanishes, 
as also clearly see in Figure [31 

We can better see what is happening by looking at 
the constraints in the Al vs Em„ plane in Figure 4. As 
we can see, there is a degeneracy between Al and Em„. 
Namely, a larger value of Y>m v decreases the lensing sig- 



nal and can be compensated with a larger Al. Since the 
SPT dataset is preferring smaller values of the lensing 
parameter, an analysis with Al — 1 forces the neutrino 
mass to be more consistent with the data. 

Is also worth mentioning that including a neutrino 
mass exacerbates the lensing problem for ACT. The lens- 
ing parameter Al is even higher when massive neutrinos 
are considered (see Table II). 



IV. DISCUSSION 

In this paper we have pointed out a tension between 
the parameter values estimated from the recent ACT 
and SPT datasets. This discrepancy, albeit not signif- 
icantly more than the 95% confidence level, is indicat- 
ing the possible presence of systematics in at least one 
of the two datasets. The SPT experiment is confirming 
the previous indications for a "dark radiation" compo- 
nent with N e s = 3.76 ± 0.34 at 68% c.l.; in particular 
we have found that N e g > 3.08 at more than 95% c.l.. 
This result is clearly interesting since, if confirmed with 
larger significance by future data, could be possibly ex- 
plained by several physical mechanisms, and would hint 
to new physics. In fact, the most conventional expla- 
nation for iVcff > 3.046 would be the presence of non- 
vanishing neutrino chemical potentials, i.e. of a cosmo- 
logical lepton asymmetry. However, as it was shown in 
Refs.@, Q through the analysis of BBN and CMB data, 
lepton asymmetries can at most account for N c g ~ 3.1, 
given the recent measurements of the neutrino mixing an- 
gle 6>i3 by the Daya Bay [H| and RENO experiments (24 1 
that exclude a zero value for #13 with high significance. 

Thus, if confirmed, a value of N g larger than 3.1 defi- 
nitely requires some non-conventional explanation. Ster- 
ile neutrinos, extra dimensions, gravity waves or non- 



6 



standard neutrino decoupling could all be viable new 
mechanisms to explain a value of N e g larger than the 
standard value (see e.g. [IS])- 

The ACT experiment is, on the contrary, fully con- 
sistent with A cff = 3.04 even when the HST and BAO 
dataset are included. There is clearly no evidence for 
dark radiation from ACT; moreover the case for a fourth, 
massless neutrino, such that N g ~ 4 is excluded at 
more than 95% c.l. when only CMB data is consid- 
ered. In particular, we found at 95% c.l. that N e g = 
2.85+°^ for WMAP9+ACT and A cff = 3.23+°^ for 
WMAP9+ACT+BAO+HST. It is interesting to notice 
that our WMAP9+ACT+BAO+HST run provides the 
constraint N c s — 3.23 ± 0.47 while a similar analysis 
from ACT gives N cS = 3.52 ± 0.39 but with Al = 1 and 
with the WMAP7 data. 

However, ACT presents a value for the lensing param- 
eter that is off by more than 95% from the expected value 
Al = 1. This result is probably more difficult to explain 
from a physical point of view than a deviation in N e f{ and 
calls for more drastic changes in the cosmological model. 
A possible way to enhance the lensing signal is to assume 
a modification to general relativity. f(R) models as those 
investigated in [6| could in principle enhance the lensing 
signal, even if it is not clear if they could enhance it by 
~ 70% and be at the same time consistent with other in- 
dependent limits coming from tests of general relativity, 
like, e.g. solar system tests. Other possible explanations 
include coupled dark energy models (see e.g. [2l| and ref- 
erences therein) . Clearly, it may be that the ACT lensing 
signal is on the contrary simply produced by some un- 
known systematic as also suggested by the inclusion of 
the ACT deflection spectrum data, that shifts the value 
to Al = 1.3 ± 0.23 ([1]). However it is not clear if this 
systematic could also affect the ACT constrain on N e s 
and other parameters. 

The SPT experiment is compatible with Al = 1 but 
is suggesting a value Al < 1 at about 68% c.l. especially 
when also the BAO and HST data are included. 

The ACT and SPT measurements of Al, even if we 
consider variation in the N c g parameter, are in disagree- 
ment at more than 95% c.l.. 

Finally, we have also considered variation in the neu- 
trino mass and show that the current indication for a 



neutrino mass from the SPT+WMAP9+BAO+HST run 
is driven by the lower lensing amplitude measured by 
SPT. If we let the lensing parameter Al to vary the indi- 
cation for a neutrino mass vanishes. Moreover, we have 
shown that the inclusion of a neutrino mass exacerbates 
the lensing problem for the ACT data with the Al even 
more discrepant with the Al — 1 case. The constraints 
on the neutrino mass from ACT are weaker when varia- 
tions in Al are considered. 

In this paper we have only considered a limited set 
of parameters but the tension between SPT and ACT 
is present also in other, relevant, parameters. The SPT 
dataset, for example, shows a preference for a negative 
running of the inflationary spectral index at more than 
95% c.l. while the ACT data is consistent with a zero 
running in between the 95% c.l. (see Figure 11 of [l[). 

We therefore conclude that the whole picture is, at the 
moment, stimulating and puzzling at the same time. The 
ACT and SPT collaborations have provided an impres- 
sive confirmation of the theoretical expectations concern- 
ing the damping tail of the CMB anisotropy spectrum. 
However they are also suggesting interesting deviations 
from the standard picture, that are unfortunately very 
different and opposite. It will be the duty of future re- 
analyses of the ACT and SPT data (possibly stemming 
from within the collaborations themselves) and experi- 
ments (e.g., Planck) to finally decide whether what ACT 
and SPT are currently seeing is due to dark radiation, 
dark gravity or more simply to an unidentified (hence, 
dark too) experimental systematic effect. 



V. ACKNOWLEDGEMENTS 

It is a pleasure to thank Andrea Marchini and 
Valentina Salvatelli for help. We also acknowledge 
CASPUR for computational support. The work of 
ML has been supported by Ministero dell'Istruzione, 
delPUniversita e della Ricerca (MIUR) through the PRIN 
grant "Galactic and extragalactic polarized microwave 
emission" (contract number PRIN 2009XZ54H2-002). 



[1] J. L. Sievers, R. A. Hlozek, M. R. Nolta, V. Acquaviva, 

G. E. Addison, P. A. R. Ade, P. Aguirre and M. Amiri 

et aL larXiv:1301.0824l [astro-ph.CO]. 
[2] Z. Hou, C. L. Reichardt, K. T. Story, B. Follin, R. Keisler, 

K. A. Aird, B. A. Benson and L. E. Bleem et al, 

larXiv:1212.6267l [astro-ph.CO]. 
[3] G. Mangano, G. Miele, S. Pastor, T. Pinto, O. Pisanti 

and P. D. Serp ico, Nucl. Phys. B 729, 221 (2005) 

[hep-ph/0506164j . 
[4] G7 Mangano, G. Miele, S. Pastor, O. Pisanti and 

S. Sarikas, Phys. Lett. B 708, 1 (2012) |arXiv: 11 10.43351 



[hep-ph]]. 

[5] E. Castorina, U. Franca, M. Lattanzi, J. Lesgourgues, 
G. Mangano, A. Melchiorri and S. Pastor, Phys. Rev. D 
86, 023517 (2012) |arXiv:1204.2510l [astro-ph.CO]]. 

[6] E. Calabrese, A. Cooray, M. Martinelli, A. Melchiorri, 
L. Pagano, A. Slosar and G. F. Smoot, Phys. Rev. D 80 
(2009) 103516 }arXiv:0908. 15851 [astro-ph.CO]]. 

[7] E. Calabrese, A. Slosar, A. Melchiorri, G. F. Smoot 
and O. Zahn, Phys. Rev. D 77 (2008) 123531 
[arXiv:0803.2309l [astro-ph]]. 



7 



Z. Hou, R. Keisler, L. Knox, M. Millea and C. Reichardt, 
larXiv:1104.2333l [astro-ph.CO] . [18] 
M. Archidiacono, E. Calabrese and A. Melchiorri, [19] 
Phys. Rev. D 84 (2011) 123008 [arXiv: 1109. 27671 [astro- 
ph.CO]]. [20] 
E. Giusarma, M. Corsi, M. Archidiacono, R. de Putter, 
A. Melchiorri, O. Mena and S. Pandolfi, Phys. Rev. D 83 
(2011) 115023 arXiv:1102.4774l [astro-ph.CO]]. 
A. G. Riess, L. Macri, S. Casertano, H. Lampeitl, 
H. C. Ferguson, A. V. Filippenko, S. W. Jha and W. Li 
et al, Astrophys . J. 730, 119 (2011) [Erratum-ibid. 732, 
129 (2011)] [arXiv: 1103.29761 [astro-ph.CO]]. 
E. Komatsu et a l, Astrophys. J. Supp. 192, 18 (2011) 
|arXiv: 1001.45381 [astro-ph.CO]] . 

K. T. Story, C. L. Reichardt, Z. Hou, R. Keisler, 

K. A. Aird, B. A. Benson, L. E. Bleem and J. E. Carl- 

strom et al, larXiv:1210. 72311 [astro-ph.CO]. 

G. Hinshaw, D. Larson, E. Komatsu, D. N. Spergel, [21] 

C. L. Bennett, J. Dunkley, M. R. Nolta and M. Halpern 

et al, larXiv: 1212.52261 [astro-ph.CO], 

A. Lewis and S. Bridle, P hys. Rev. D 66, 103511 (2002) [22] 
[arXiv:astro-ph/0205436] . 

Beutler, F., et al., Montly Notices of the Royal Astro- 
nomical Society, 416, 3017, 2011 [23] 
Padmanabhan, N., Xu, X., Eisenstein, D. J., Scalzo, R., 
Cuesta, A. J., Mehta, K. T., & Kazin, E. 2012, ArXiv [24] 



e- prints, arXiv: 1202.0090 

L. Anderson, et al. 2012, ArXiv e-prints. larXiv: 1203.65"MI 
C. Blake et al., Montly Notices of the Royal Astronomical 
Society, 425, 405, 2012 

S. Hannestad, A. Mirizzi, G. G. Raffelt and 
Y. Y. Y. Wong, JCAP 1008 (2010) 001 |arXiv:1004.0"695l 
[astro-ph.CO]]; A. Melchiorri, O. Mena and A. Slosar, 
Phys. Rev. D 76 (2007) 041303 |arXiv:0705.2695l 
[astro-ph]]; T. L. Smith, E. Pierpaoli and 
M. Kamionkowski, Phys . Rev. Lett. 97 (2006) 021301 
[arXiv:astro-ph/0603144] ; P. Binetruy, C. Deffayet, 
U. Ellwanger and D. Langlois, Phys. Lett. B 477 
(2000) 285 [arXiv:hep-th/9910219l ; T. Shiromizu, 
K. i. Mae da and M. Sasaki Phys. Rev. D 62 (2000) 
024012 |arXiv:gr-qc/9910076j ; S. Mukohyama. Phys. 
Lett. B 473, 241-245 (2000), doi = "10.1016/S0370- 
2693(99)01505-1" 
M. Martinelli, 
and O. Mena, 



L. Lopez Honorez, 
Phys. Rev. D 81 



A. Melchiorri 
(2010) 103534 



arXiv:1004.2410 [astro-ph.CO]]. 
J. Dunkley, E. Calabrese, J. Sievers, G. E. Addison, 
N. Battaglia, E. S. Battistelli, J. R. Bond and S. Das 
et al.. larXiv:1301.0776l [astro-ph.CO], 
F. P. An et al. [DAYA-BAY Collaboration], Phys. Rev. 
Lett. 108, 171803 (2012) |arXiv:1203. 16691 [hep-ex]]. 
J. K. Ahn et al. [RENO Collaboration], Phys. Rev. Lett. 
108, 191802 (2012) |arXiv: 1204.06261 [hep-ex]]. 



